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Abstract 

The role of thermal fluctuations on the conformational dynamics of a single closed filament is 
studied. It is shown that, due to the interaction between charges and bending degrees of freedom, 
initially circular aggregates may undergo transformation to polygonal shape. The transition occurs 
both in the case of hardening and softening charge-bending interaction. In the former case the 
charge and curvature are smoothly distributed along the chain while in the latter spontaneous kink 
formation is initiated. The transition to a non-circular conformation is analogous to the phase 
transition of the second kind. 

PACS numbers: 87.15.-v, 63.20.Pw, 63.20.Ry 
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I. INTRODUCTION 



Conformational flexibility is a fundamental property of biological systems which deter- 
mines their functioning; Even modest conformational changes modify long-range 
electronic interactions in oligopeptides 0] , they may remove steric hindrances and open the 
pathways for molecular motions which are not available in rigid proteins [5|. The DNA 
conformation in the nucleosome core is crucial for gene replication, transcription and re- 
combination [(J. Recent DNA cyclization experiments have shown the facile in 
vitro formation of DNA circles shorter than 30 nm (100 base pairs) which is even shorter 
than commonly accepted persistence length 50nm (150 base pairs). This means that the 
worm-like chain model does not work for such short DNA molecules and to explain this 
phenomenon one should allow local softenings of DNA which facilitates disruptions (kinks) 
in the regular DNA structure 1(3, ll|, [3]. According to [lo| the kink formation is due to 
strong DNA bending while in [9] that the softening originates from Watson-Crick base-pair 



breathing. An alternative approach which allows to avoid kinking was proposed in Ref. 13] 
where a class of models with nonlinear DNA elasticity was introduced. It was shown in 13] 
that a "subelastic chain" model, in the frame of which the bending energy is proportional 
to absolute value of curvature, can reproduce the main features of Cloutier and Widom's 
experiments jj]. 

Quite recently a simple, generic model for charge-curvature interactions on closed molec- 
ular aggregates was proposed [jjl]. It was shown that the presence of charge modifies (softens 
or hardens) the local chain stiffness. It was found that due to the interaction between charge 
carriers and the bending degrees of freedom the circular shape of the aggregate may become 
unstable and the aggregate takes the shape of an ellipse or, in general, of a polygon. It 
was shown also that when the charge- curvature interaction leads to softening the local chain 
stiffness kinks spontaneously appear in the chain. 

These results were obtained by using the mean-field approach where thermal fluctuations 
are ignored and strictly speaking, this approach is valid only for zero-temperature. In the 
case of finite temperature the interaction with environment and thermal fluctuations have 
to be considered. 

The aim of this paper is to extend the results of Ref. [yj to the case of finite temperature. 
We study the charge- induced conformational transformations of closed molecular aggregates 
in the presence of thermal fluctuations which we model in the frame of Langevin dynamics. 
The paper is organized as follows. In Sec. II we describe a model. In Sec. Ill we present 
an analytical approach to the problem. In Sec. IV we display the results of numerical 
simulations and compare with the analytical results. In Sec. V we discuss some concluding 
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remarks. 



II. THE MODEL 

We consider a polymer chain consisting of L units ( for DNA each unit is a base pair) 
labelled by an index I, and located at the points r*/ = (xi, yi), I = 1 . . . L. We are interested 
in the case when the chain is closed, therefore we impose the periodicity condition on the 
coordinates, 

n = n +L . (l) 

We assume that there is a small amount of mobile carriers (electrons, holes in the case of 
DNA, protons in the case of hydrogen bonded systems) on the chain. The Hamiltonian of 
the system can be presented as the sum 

H t ot — H + H stoch . (2) 

The first term in this equation is the Hamiltonian of an isolated filament introduced in Ref. 

H 

H = U b + U s + H el + H e i^ con f. (3) 

Here 

U b = ~y2 K[ (4) 

is the bending energy of the chain where 

«j = \ti -ti-i\ = 2 sm- (5) 

determines the curvature of the chain at the point /. Here 

ri+i - fi 

Tl = p; ZTT 6 

\n+i - ri\ 

is the tangent vector at the point I of the chain and a\ is the angle between the tangent 
vectors fj and r/_i, k is the elastic modulus of the bending rigidity (spring constant) of the 
chain. The term k^/k^^ in Eq. (Tj0) gives the penalty for too large bending deformations. 
Here the parameter n max = 2 sin (a max /2) is the maximum local curvature with a max being 
the maximum bending angle. The second term in Eq. (J3j) 

i 

determines the stretching energy with a being an elastic modulus of the stretching rigidity 
of the chain and a is the equilibrium distance between units (in what follows we assume 
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a = 1). We take the simplest theoretical model for charge carriers, a nearest neighbor tight 
binding Hamiltonian in the form 

H el = jJ2 [Pi - VHi , (8) 
i 

where ipi is the wave function of carrier localized at r\ and J measures the carrier hopping 
between adjacent sites. The last term in Eq. Q represents the charge-curvature interaction. 
In the small curvature limit it has the form 

H e l- con f = -- J^X l^l 2 + «J-l) , (9) 

I 

here \ is the coupling constant. Combining Eqs. (Hj) and (J9j), we notice that the effective 
bending rigidity changes close to the points where the electron (hole) is localized. For 
positive values of the coupling constant \ there is a local softening of the chain, while for x 
negative there is a local hardening of the chain. 
The quantity 



-yEl^l 2 (10) 



i 

gives the total density of charge carriers which can move along the chain and participate 
in the formation of the conformational state of the system. The second term in Eq. §2§ 
describes the interaction of the filament with fluctuating environment 

H stoch = J2rrRi(t) (11) 

i 

where the stochastic forces Ri(t) = (Xi(t), Yi(t)) are the Gaussian white noise 

(x t (t)) = m)) = o, 

= (Y l (t)Y l/ (t>)) = 2D5 ll ,5(t-t'), (12) 
(X l (t)Y v (f)) = 

with D the standard deviation. 

To analyze the evolution of the shape of the filament, it is convenient to introduce the 
radius-of-gyration tensor I as in Ref. HQ. Its components are 

I M (t) = ^^(xKt)-z c (t)) 2 , 
i 

W*) = ^E(w(*)-y c (*)) 2 , 
i 

K(t) = jJ2 ( X M ~ * C W) " y C H)) (13) 
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where 



x 



Xt),y c (t)) = jJ2(xi(t)Mt)) 



(14) 



is the center-of-mass coordinate. The square roots of the two eigenvalues R q = y I g , q = 1,2 
of the tensor I give the two principal radii of the system. They express the sizes of the 
filament along the major and minor axis. As it is seen from Eqs. ffl3|) the eigenvalues have 
the form 

' 1 (15) 



L l,2 



I, 



I 



+ 4I 2 

1 xy 



Indexes 1 and 2 correspond to the + and — sign respectively. To characterize the shape of 
conformation it is convenient to introduce the quantity 



A = h-h 



I - I ) 2 + 41 



xy ■ 



(16) 



defined as the "aspherity" 16]. It characterizes the shape's overall deviation from circular 



symmetry which corresponds to A = 0. 



III. ANALYTIC APPROACH 

The aim of this section is to develop an analytical approach, which provides a better in- 
sight into the physical mechanism of conformational transformations induced by the charge- 
curvature interaction in the fluctuating media. We will assume that the characteristic size of 
the excitation is much larger than the lattice spacing and replace ipi and by the functions 
ip(s,t) and f(s,t), respectively. Here the arclength s is the continuum analogue of I. We 
assume that the chain is inextensible and this assumption is expressed by the constraint 

\d s ^\ 2 = 1. (17) 

which is automatically taken into account by choosing the parametrization 

d s x(s) = sin6>(s), d s y(s) = cos9(s) (18) 

where the angle 9(s) satisfies the periodicity condition of Eq. ([Tj) 

6(s + L) = 2tt + #(s), (19) 

L L 

J cos 6{s) ds = J sin 6(s) ds = 0. (20) 
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In the frame of the parametrization (1TB]) the shape of the chain is determined by the equations 

s s 

x(s) = J sin 9{s') ds', y(s) = J cos 9{s') ds' . (21) 



In the continuum limit the curvature ([5]) takes the form k(s) = \dgf(s)\, which is given by 

k(s) = d s 9. (22) 
The continuum version of the total Hamiltonian of the system can be written as the sum 

H t ot = H + H stoch {t) (23) 

where 

H = j\j\dM 2 + (|-xhA 2 |) (d s 9) 2 \ds (24) 
o ^ ' 
is the analogue to the Hamiltonian ([3]) and 

L 

Hstoch(t) = J {X{s, t) x{s, t) + Y(s, t) y(s, t)) ds (25) 



gives the interaction of the chain with the fluctuating environment (jTTj) in the continuum 
limit. The stochastic forces R(s,t) = (X(s,t), Y(s,t)) are the continuum version of the 
forces Xi(t) and Yi(t). They obey the relations 

(X(s,t)X(s',t')) = (Y(s,t)Y(s',t')) = 
2D5(s-s')8(t-t'), 

(X{s,t)Y(s',t'))=0. (26) 

We will restrict our analysis to the case when the filament shape only slightly deviates from 
the circle. Therefore in the derivation of the Hamiltonian ([53]) we could neglect the term 
K 2 /K>max m ^e denominator of Eq. (jlj). 
By using the Madelung transformation 

V(M) = v'pM^''" (27) 

where p(s,t) is the charge density and (fi(s,t) is the phase, the Hamiltonian (1241) can be 
written as follows 



H 



L 

/] 









Ap 



2 



P (<9 S 0) 



\-xp) {d s ef\ds. (28) 
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The dynamics of the system is governed by the Hamilton equations for the charge variables 
p(s,t), <f)(s,t) 



6 -f = and ^=0 (29) 
dp 8<p 



where 



C 



pdt4>ds — H 



(30) 



is the Lagrangian of the system and 8/8(-) is a variational derivative. By introducing the 
dissipation function 



T = r) i J {d t rf ds 



(31) 



the Langevin equation for the position f(s, t) can also be written in the variational form 

8 J- 8 Htot / qo n 

8d t f{s,t) = ~8f{s,t)' ^ ' 

Eqs. ( ]29~|) and fl32|) should be considered with the periodicity condition (IT91 and the closure 
condition (1201 for the position r(s, t). Now the periodicity conditions for the charge variables 
take the form 

p(s)=p(s + L), cf>(s)=<f)(s + L) (33) 
and the normalization condition in ffTOl becomes 



pds = v. 



(34) 



To take into account the periodicity conditions (|T9|) and (|33|) we can expand the curvature 
and the charge variables in the Fourier series: 



p{s,t) = V 



2 7TJ S 



l+E&Wcos(^ 

4^0 \ 



S,t) = COS f 

i>2 ^ 



2 7TJ S 



(35) 

(36) 
(37) 



Note that the first harmonic with j = 1 does not contribute to the Fourier expansion 
fl35|) due to the closure condition (TSUI) . The coefficient z/ in the expansion (136]) takes into 
account the normalization condition (|34l) . Different harmonics in the Fourier expansion (|35|) 
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represent different types of shape deformation. For example the term with j = 2 determines 
an elliptic deformation, the term with j = 3 represents a triagonal deformation, etc (see 
(l^ for more detail). For the sake of simplicity we consider only the elliptic deformation of 
the filament (j = 2) and restrict ourselves to the case when the deviations from the circular 



shape are small and the charge distribution along the chain is smooth: \c2\1 \p2\1 



-3\i 



I Pi I 



<C 1 for j > 3. Thus the expansions (13"5]) - (I3"T|) reduce to 



d s e(s,t) 



2tt 



p(s,t) = v 



1 + c{t) COS 

1 + ((t) COS 



An s 
Air s 



(s,t) = $(t) cos 



L 

An s 



A < 

(38) 
(39) 
(40) 



where we omitted subscript in the notations for the Fourier harmonics. 
Inserting Eqs. (1581-fllOl into Eqs. (125) . (125]), (1501) and (151]) we get 



eff 



L d§ 



n. 



Here 



H 



2 

^- (-4Ji/y/l-£ 2 + 8 Jz/$ 2 + k e c 2 -A X v(> 



(41) 



(42) 



is the effective Hamiltonian with some irrelevant constants being omitted. In Eq. (f42|) 
k e = k — 2 x v is an effective bending rigidity of the filament for the case when the charge is 
unformly distributed along the chain. 

J? 

7~(-stoch{t) 



An 2 



+Y 



L_s 
2^T 





s 

, ( )/cos( S ' + |sin(2 s0 )^ 



(43) 



ds 



is the effective interaction with stochastic forces, and 

"-'Mis 

is the effective dissipative function. The damping coefficient 6(c) has the form 

2tt 



(44) 



b(c) 



1 L 3 

I" 7 (2tt) j 



(o; 2 (c, s) + a 2 (— c, s)) (is 



(45) 
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where the notation 



a{c, s 



s 

/ sin (2 s') sin (V + | sin(2 s')) ds 



(46) 



is introduced. Note also that in the derivation of Eq. (I4"51) we took into account the period 
icity in the stochastic terms, X(s + L,t) = X(s, t) and Y(s + L,t) = Y(s, t). 

Equations of motion for the quantities $, ( and c follow from Eqs. (l28l) . (J29l) . ( l4Tj) - fl44l 
and they have the form 



rfC 32tt 2 



dt L 2 



dc 



2vr 2 



dt 6(c) L 



where 



2tt 



/(c,t) 



46(C) 7T 2 



V ( ) «(-c,^i 





- V I ^,t) n(r..s) 

2 7T 



(47) 

(48) 
(49) 



(50) 



is an effective stochastic force. 

Note that, in terms of the Ansatz (138j) . the aspherity A, which is defined by Eq. (fl6l) 
can be written approximately as 

/- 

(51) 



A = — c. 



8tt 2 

Let us analyze the cases of zero temperature and finite temperature separately. 



A. Deterministic behavior: zero temperature limit 

In the no-noise case the dynamics of the system is described by Eqs. (j4Tj)-(l4"9l) with 
X = Y = 0. The system under consideration is characterized by the control parameter 

e = ^ (52) 

which is the ratio of the deformation energy (i.e. the energy shift due to the charge-bending 
interaction) 

2 2 

XT 



^ (53) 
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with respect to the dispersion energy 

&dis P = \ (jpj vJ. (54) 

A simple analysis shows that when the charge-curvature coupling is weak such that the 
control parameter, £ < 1, these equations have a unique stationary point $ = 0, ( = and 
c = 0. This state corresponds to a uniformly distributed charge along the circular filament. 
When £ > 1 there are two equivalent stationary states 

$=o ' c±=± \A^' c±= ^~ c± (55) 

which represent an elliptically deformed filament with a spatially non-uniformly distributed 
charge. The two solutions c± correspond to two mutually orthogonal directions in which the 
filament may be elongated. Note that in the case of softening charge-curvature interaction 
(x > 0) the maxima of the curvature and the charge density coincide while in filaments with 
hardening charge-curvature interaction (x < 0) the curvature of the filament is minimal (the 
filament is locally more flat) in the places where the charge density is maximal. 



B. Charge-charge correlation effects 

The aim of this subsection is to clarify the role of interaction between charge carriers in 
the formation of polygonally shaped aggregates in the zero-temperature limit. In describing 
the charge-charge repulsion effects, we will use a on-site interaction in the form 

H el _ el = l -VY, M (56) 

n 

which in the continuum limit in terms of charge variables ( J27l) reads 

L 

Eei-ei = \v J P 2 ds. (57) 



The parameter V in Eqs. (156]) and fl571) characterizes the strength of the interaction. Thus 
the Hamiltonian of the system with account of charge-charge interaction effects has the form 

L 



H CC = H + -V I p 2 ds. (58) 



o 
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where the Hamiltonian H is given by Eq. (1281) . Inserting in Eq. (1581) the Ansatz (|38|) . we 
get 

8 Jz/$ 2 + /c e c 2 -4x^Ccj. (59) 

An inspection of the function (1591) shows that in the zero-temperature limit the spatially 
uniform charge distribution along the circular filament becomes unstable for the control 
parameter (l52i) satisfying the inequality 

L 2 V 

(> l + (60) 

Thus the elliptic shape more easily arises in short filaments with strong charge-bending 
interaction x an d relatively weak interaction between charge carriers V. 

In what follows we will assume that the charge-bending interaction is strong, A^ef 3> 
v 2 V 1 4, and for the sake of simplicity will neglect the interaction between charges. 

C. Stochastic behavior: finite temperature 

We will study the role of thermal fluctuations by using a formalism of the Fokker-Planck 
equation. To this end we introduce the probability distribution density 

P (c, C, t) = (s{c- c{t)) 5 (C - C(t)) S ($ - $(t)) ) (61) 

As it is seen from Eqs. ( |26l) and ( |50l) the stochastic forces f(c,t) represent Gaussian white 
noise with the mean value 

(f(c,t)) = (62) 

and the two time covariance given by 

(f(c,t)f(c> } t')) = 2F(c,c')6(t-t') (63) 

where 

3 



_ d l a r 

(C ' C) ~ 6(c) b(d) 8^ J 





+a(c, s) a(c', s) ) ds. (64) 
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It is straightforward to obtain (see e.g. jl7[]) that the Fokker-Planck equation which describes 
the time evolution of the probability distribution ( loll) of a set of Langevin equations f j47l) -( )49|) 
in the Stratonovich sense has the form 

P 



d t P = -d c (P d$H) + 3$ (P d c H) + d c 



b(c 



P (d c F(c,c')) 



d 2 c (PF(c,c)) 



where the Hamiltonian TC is given by Eq. fj42l Inserting into Eq. fl65l) the relations 



F(c,c) 



T 



, d c F(c,c') 



T d 



2 dc \b(c 



1 



(65) 



(66) 



which follow from Eqs. (J4"5j) . fj46l) and ([M]) . we obtain the equation for the probablity 
distribution P (c, (, $; t) in the form 

P 



d t P = -d c (P d^rt) + <9$ (P 9 C H) + d, 

1 





2 



pi 



(ic \6(c) 



b(c) 



b(c 



P 



(67) 



It is interesting to note that if the stochastic force f(c,t) in the Langevin equations 
(HZD-dUD is replaced by 



mod 



where £(i) is a white noise with 



(68) 



(69) 



We obtain the same Fokker-Planck relation for this new set of Langevin equations. This 
means that the stochastic equations P71) - (l49l) choosing f(c,t) from Eqs (12"6"1) and (|50|) or 
from Eqs. fl68|) and fl69l) are equivalent. 

The stationary probability distribution (i.e. the solution of the Fokker-Planck equation 
(J6"7|) for t — > oo) is given by 



P 



c e 



(70) 



where C is a normalization constannt. 

Assuming that c < 1, one can expand the function a(c,s) from Eq. (j46p into a Taylor 
series and keeping only leading terms one can obtain that 



In 6(c) w const + 0.04 c 2 . 



(71) 



Combining Eqs. fj42l) . flTOj) and fl7T|) we obtain that the effective bending rigidity (i.e. the 
coefficient in front of c 2 in the expression Tl+ 1 In 6(c)) becomes temperature dependent and 
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takes the form k e — 0.002L T. We will assume that the persistence length l p = k e /T satisfies 
the inequality l p ^> 0.002L and in what follows we will neglect the temperature dependence 
of the effective bending rigidity. 

By integrating the stationary probability distribution (1701) over the charge variables $ 
and ( (note that £ 6 (—1,1) because under definition the charge density p(s,t) > 0) we 
obtain a reduced distribution 



tt/2 

V(c) = J dO sin 6 e 13 sinf> cosh (j3^ c cos ej 



exp 
where 



2£ Vj' 



(72) 



tt/2 

AT 



(73) 



cxp 



d c J d9 sin sin9 cosh c cos ^ 

2^ vjy 

is the normalization constant. In Eqs. (1721) and (1731) /3 = 4 7r 2 J v/(LT) is a dimensionless 
inverse temperature. The function V(c) gives the probability of finding the curvature c in 
the interval (c,c + d c) irrespective of the magnitude of the charge variables. 

There are two areas in the parameter space £) where the probability density (1721) as a 
function of the curvature parameter c behaves qualitatively different as is shown in Fig. HJ 
These two areas are separated by the curve 

t = P 2 +7T (/!(/?) 

^ vr (7 (/3) + L (/3)) - 2 (JiOS) + LxCS)) 1 J 

where I n ((3) is the Bessel function of imaginary argument and L n (/3) is the Struve function 
jisl ]. Below this curve (ie. in the unshaded area of the phase diagram presented in Fig. [1]) 
the distribution is single-modal. The probability density ( 1721) in this case has a maximum 
at c = (see Fig. [2] when T = 0.3). This means that the most probable conformation 
state of the filament is a circle. Above the curve (1741) (i.e. in the shaded area of Fig. 
H]) the probability distribution is bimodal. The function ( 1721) has two equivalent maxima 
±c m (see Figf2] when T = 0.05 and T = 0.1). In this case the most probable state is 
an elliptically deformed filament. Two maxima correspond to two mutually orthogonal 
directions of elongation. 

It is straightforward to obtain that in the vicinity of the curve given by Eq. (1741) the most 



probable value of the chain curvature is determined by the expression c m = 7(£) y/T c — T 
where T c (£) is a critical temperature (i.e. the solution of Eq. (174|) ) and 7(£) is some 
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1.5 



. 5 




0.1 0.2 0.3 0.4 0.5 (5 

FIG. 1: The probability distribution phase diagram. In the shaded area the probability distribution 
is bi- modular and in the unshaded area it is unimodular. 
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FIG. 2: The probability density distribution of c variable for J = k e = 1, v = 0.25, x = 2, o = 10 5 
and three different values of the temperature: T = 0.05 (dashed curve), T = 0.1 (thin curve), 
T = 0.3 (thick curve) 



coefficient. Thus, in the framework of our Ansatz (|38|) the transition of the filament in the 
fluctuating environment to a non-circular conformation may be considered as a noise- induced 
phase transition of the second kind. 

By using the probability distribution (1721) and the expression (l5Tl) for the filament as- 
pherity, one can calculate the equilibrium value of the normalized aspherity 

(A) Mx(T) 



(Ao) M^O)' 

and the equilibrium value of the relative standard deviation of the aspherity 



(75) 



(A) 



M 2 (T) 
Ml(T) 



(76) 
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where M n (T) are the moments defined as 

1 

M n (T) = J V (c) \c\ n dc, n = l,2,... (77) 



IV. NUMERICAL STUDIES 

The dynamics of the filament is described by the Schrodinger equations 

^ = -^w (78) 

and the Langevin equations 

d dH -> , . . . 

"s f ' = -^ +fliW - (79) 

with the Hamiltonian H being defined by Eq. ([3]). Thus the conformational dynamics is 
considered in an overdamped regime with r\ being the friction coefficient. In accordance 
with the fluctuation-dissipation theorem the standard deviation D is proportional to the 
temperature T, thus D = rj T. 

The set of stochastic differential Eqs. (ITS]) . (1791) is solved numerically by the use of an 
implicit Euler method with a = 0.5 as implicitness parameter, which is the trapezoidal 
rule. The integration of the stochastic term is done by using the strong Taylor scheme of 



first order described in Ref. [19|. The time step chosen for running the simulations was 
At = 10~ 2 . To verify the precision of the results we compared with data obtained for 
different time steps. The value of the position and charge of Z'th particle at time t n = nAt is 
denoted Z l n = (xi, yi, (fi) n , F(Z n ) denotes the deterministic part of Eqs. (J7Hj) . (J79]) and W n = 
(Xi(t n ),Yi(t n ),0) is the corresponding component of the white noise. Thus the numerical 
scheme becomes 

Z l n+1 = Z l n + [aF(Z n+1 ) + (1 - a)F(Z n )]At 

+W n VAi. (80) 

The above system of nonlinear equations is implicit and it is solved by a hybrid method 
provided by the minpack FORTRAN library and the random numbers present in the white 
noise are generated by the ranlib library, both accesible from the "netlib" repository in Ref 



20]. 



In our paper we were mostly concerned with the role of thermal fluctuations in the 
process of the shape transformations. Therefore in our simulations we have chosen a set of 
parameters for which in the zero-temperature limit the most energetically favorable state 
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is an elliptically deformed filament. Without loss of generality, in this section we show the 
results of the numerical simulations produced for L = 36 units and charge density v = 0.5. 
A system of this size provides a clear visualization of the properties of the model and does 
not demand too heavy computational time. The systems of the same size but with smaller 
values of the charge density (however, still inside the area where the elliptic state is stable) 
in the presence of thermal fluctuations require more time to reach an equilibrium state. 
As initial condition for electric charge density ipi we use the same magnitude at all points, 
corresponding to an equally distributed charge density. Initially, all the lattice points were 
symmetrically distributed along the circle of an appropriate radius. In what follows we have 
chosen the damping coefficient r\ and the bending rigidity k equal to unity. Moderate changes 
in these parameters do not modify significantly the dynamics of the system. To avoid big 
stretching of the nearest bonds of the chain we fix the parameter a = 10 5 . We considered 
both the cases of the hardening and of the softening electron-curvature interaction. 

A. The case of hardening charge- curvature interaction 

Typical final shapes of the filament for three different values of the noise intensity D = 
0.01, D = 0.1 and D = 0.2 are shown in Fig. [3]where the averaged aspherity decreases as the 
noise increases. Fig. H] shows the time evolution of the aspherity A for three different values 
of the noise D. The left panel shows the overall behavior of this quantity (including transient 
processes) while the right panel presents its steady state evolution. As it is seen from Fig. 
H]the thermal fluctuations excite the dynamics of the system and facilitate the transition to 
an anisotropic state. Clearly the transient period shortens when the temperature increases. 

It is remarkable that in the case of weak noise, D = 0.01, in the aspherity time evolution 
there exists a small plateau for 600 < t < 1000 (see Fig. Hj). This plateau corresponds to an 
intermediate conformation of the filament when it takes a triangular shape (see Fig. [5]). 

The mean value of the saturation aspherity (A) after transient time ti which we define 

as 



ti 

decreases as the temperature increases (see Fig. [6]) and the filament takes on a less 
anisotropic shape while the relative standard deviation 




(81) 




(82) 
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FIG. 3: The top panel shows the equilibrium shape of the chain (solid line) against the initial 
circular shape (dashed line). The bottom panel shows the charge distribution \t(j\ 2 (solid line) and 
curvature k (dashed line) for different noise intensity: D = 0.01 (a), D = 0.2 (b), and D = 0.4 
(c) in the case of hardening, v = 0.5, \ = ~~ 4, <r = 10 5 and J = 0.25 at time t = 4000. The grey 
shadow represents the local charge density in the chain. 
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FIG. 4: Aspherity versus time for different noise intensity, D = (thin line), D = 0.05 (thick line) 
and D = 0.25 (dotted line). The parameters used are v = 0.5, x = — 4, a = 10 5 and J = 0.25. (a) 
Full time simulation < t < 4000, (b) Detailed behavior after transient time 2000 < t < 4000 



is an increasing function of temperature. Here the aspherity variance is given by the expres- 
sion 



t 2 

((AAf) = ^-L- J (A(t) - (A)) 2 dt. 



(83) 
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FIG. 5: The intermediate shape of the filament at time t = 500 (left panel) and the corresponding 
charge distribution \ip\ 2 (solid line) and curvature k (dashed line) along the chain (right panel) for 
D = 0.01. The parameters used are the same as in Fig. [3] 

The shape of the aggregate is well-defined when yj ((AA) 2 ) / (A) <C 1. There is a critical 
value of the noise intensity D cr when the relative standard deviation reaches the value 1/2. 
Then for D > D cr the shape fluctuations are so strong that they become chaotic. For 
v = 0.5 this value is approximately D cr rs 0.4. 

Thus one can conclude that in the case of hardening charge-bending interaction the 
mean-field approach introduced in [14j works well in the weak- noise limit: 

(i) The ellipse-like conformation is the equilibrium state of the shape evolution. 

(ii) In the equilibrium state the charge distribution is nonuniform along the chain and 
the charge is concentrated in the places where the filament is more flat. 

(iii) The triangular conformation of the filament exists but it is a metastable state. 
Strong noise (i.e. high temperature) changes qualitatively the conformational dynamics 

keeping the shape of the filament almost circular, while in the zero-temperature limit it 
takes on an elliptic shape. 

The corresponding quantities are shown in Fig. [6j Comparing these relations with the 
numerical data, also presented in this figure, one can conclude that there is a qualitative 
agreement between the analytical results obtained in the framework of the simple Ansatz 
( I38p and the results of numerical simulations. 



B. The case of softening charge- curvature interaction 

By carrying out the simulations described for the hardening charge-bending interaction 
(x < 0) we took into account that the presence of charge increases the bending rigidity 
locally and the large bending deformations do not occur. Therefore there was no necessity 
to introduce any bending energy penalty and in Eq. (j3J we put n max — > oo. However, 
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FIG. 6: The normalized mean value of the aspherity (A) / (Aq) (thick solid line) and its relative 
standard deviation \l ((Avl) 2 ) /(A) (thick dashed line) versus noise from the simulations and as 
stationary solution of the Fokker-Planck equation (thin solid line and dotted line respectively) in 
the case of hardening, v = 0.5, x = — 4, a = 10 5 and J = 0.25. 

the filaments with softening charge-bending interaction are much more flexible and large 
deformations with the bending angle a.\ > tt/2 appear rather easily. To avoid excessive 
bending we carried out our numerical simulations under the assumption that the maximum 
possible bending angle a max is equal to ir/2. This means that we let n max = \/2 in Eq. (j4j). 
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FIG. 7: Same as Fig. Hfor the case of softening, v = 0.29, x = 2, a = 10 5 and J = 0.22. (a) Full 
time simulation < t < 4000, (b) Detailed behavior after transient time 3000 < t < 4000 

Comparing the time evolution and the mean value of the aspherity in the case of softening 
charge-bending interaction presented in Fig. [7] and Fig. [8] with the case of hardening charge- 
bending interaction (Fig. [Hand Fig. [6]) we see that they are qualitatively similar. However, 
the equilibrium shapes of the filaments (Fig. [9]) and the charge distribution (Fig. [TUI) after 
transient time t = 4000 differ drastically. The softening charge-bending interaction initiates 
kink formation while smooth shapes of filaments are characteristic for the case of hardening 
charge-bending interaction. Almost all the charge is concentrated now in the areas of the 
kinks. It is interesting to notice that small and moderate thermal fluctuations facilitate the 
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FIG. 8: The normalized mean value of the aspherity (^4) / {Aq) (solid line) and the relative standard 
deviation \J ((AA) 2 } / (A) (dashed line) versus noise for the case of softening with v = 0.29, % = 2, 
a = 10 5 and J = 0.22. 

kink formation while for sufficiently high magnitudes of the noise intensity filaments takes 
on almost circular shape and the charge distribution is more uniform along the filament (see 
Fig. E^d) and Fig. IPOlfd)). 

y 




y y 




X 



FIG. 9: Equilibrium configurations of the chain in the case of softening, v = 0.29, x = 2, cr = 10 5 
and J = 0.22, for noise intensities: (a) D = 0.05, (b) D = 0.10, (c) D = 0.20 and (d) D = 0.30 
(solid line) against the initial circular conformation (dashed line). The grey shadow represents the 
charge density. 

We carried out a few runs of numerical simulations for a set of parameters which satisfy 
the inequality (!60l) . We found out that for small enough V the filament shape behavior and 
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FIG. 10: Equilibrium charge distribution \ip\ 2 (solid line) and curvature k (dashed line) along the 
chain for the same parameters as in Fig. [9j 

aspherity behavior are essentially the same as in the case of no-interaction between charges. 
More detailed studies of effects of charge-charge interactions will be published elsewhere. 



V. CONCLUSIONS AND DISCUSSION 

In this paper, we have investigated the role of thermal fluctuations on the charge-induced 
conformational transformations of closed semiflexible molecular chains. We have found that 
the results obtained in the mean-field approach [3] are rather robust in the systems where 
the presence of charge hardens the local chain stiffness, the charge-curvature interaction 
counteracts the collapse of the chain and the mean-field picture survives. In the presence of 
white noise when the charge density and/or the strength of the charge-curvature coupling 
exceed a threshold value, the spatially uniform distribution of the charge along the chain and 
the circular, cylindrically symmetric shape of the chain become unstable. In this case the 
equilibrium state of the system is characterized by a spatially nonuniform charge distribution 
along the chain which takes on an ellipse-like form. The transition to an anisotropic spatially 
nonuniform conformation is analogous to the phase transition of the second kind in the 
condensed matter physics. 

In the case of hardening charge-bending interaction the charge and curvature distribution 
along the filament are smooth while for softening charge-bending interaction there sponta- 
neously create kinks where the smooth filament structure is disrupted. Almost total amount 
of excess charge is concentrated in the vicinity of the kinks. 
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